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ABSTRACT 


Experimental  heat  transfer  studies  have  been  performed  on  three 

screens  simulating  a  porous  parachute.  The  average  Nusselt  number  was 

measured  as  a  function  of  the  sonic  Reynolds  number.  The  range  of  the 

sonic  Reynolds  number  extended  from  5  to  about  1400.  Three  different 

porosities  and  two  different  pressure  ratios  were  used.  Within  the 

range  of  porosities  and  pressure  ratios  investigated,  no  significant 

effect  of  these  parameters  on  the  average  ne3t  transfer  coefficient 

could  be  established.  The  Nusselt  numbers  for  screens  are  shown  to  have 

a  behavior  quite  similar  to  the  Nusselt  number  for  a  single  cylinder  in 

x 

cross  flow.  The  relation  Nu(L  +  D)/Re*  =  0.5  appears  to  be  an  acceptable 
representation  of  the  average  Nusselt  numbers  on  screens  for  sonic  Rey¬ 
nolds  numbers  above  400  for  which  rarefaction  effects  are  not  present. 

The  distribution  of  this  Abstract  is  unlimited. 


iii 

(reverse  is  blank) 


TABLE  OK  CONTENTS 


Section  Page 

I  INTRODUCTION  1 

IT  r utuncmAL  KAUlMlt  6 

A.  Low  Density  Wind  Tunnel  6 

B.  Test  Section  16 

C.  Instrumentation  19 

III  EXPERIMENTAL  STUDIES  21 

IV  DISCUSSION  OF  RESULTS  29 

V  CONCLUSIONS  39 

APPENDIX  I.  DATA  REDUCTION  TECHNIQUE  40 

APPENDIX  II.  ESTIMATE  OF  CONDUCTION  LOSSES  42 

REFERENCES  44 


ILLUSTRATIONS 


Figure 

1.  Sketch  of  a  Porous  Parachute 

2.  Photograph  of  the  Wind  Tunnel 

3.  Sketch  of  Experimental  Facility 

4.  Exploded  View  of  Test  Chamber 

5.  Aluminum  Plate  with  Test  Section  and  Housing  Blocks 

6.  Aluminum  Plate  with  Nozzle 

7.  One  Test  Screen 

\ 

8.  Example  of  Extrapolation  Technique  Used  to  Obtain  R^ 

9.  P  /P,  versus  Re*  for  Screen  1 

up  down 

10.  P  /P,  versus  Re*  for  Screen  2 

up  down 

11.  P  /P.  versus  Re*  for  Screen  3 

up  down 

12.  Nusselt  number  versus  Re*  for  Screen  1 

13.  Nusselt  number  versus  Re*  for  Screen  2 

14.  Nusselt  number  versus  Re*  for  Screen  3 

15.  Comparison  of  Present  Results  with  Results  for  a  Cylinder 

,  ,  1/2 

16.  Nu/ (Re*)  versus  Re*  for  All  Three  Screens 

17.  Comparison  of  Experimental  Results  with  Theoretical  Prediction 

18.  Comparison  of  Experimental  Results  and  the  Results  Corrected 
for  Estimated  Heat  Conduction  Losses 


Page 


2 

7 

8 


13 

14 
17 
24 
26 

27 

28 

30 

31 

32 


33 


37 


39 


43 


vi 


Reproduced  Prom 
Best  Available  Copy 


SYMBOLS 


d 
li 
I 
k 

K.n 
Ma 
Nu 

r 

up 
^down 


Po 

q 

R 
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Superscripts 


*  denotes  sonic  conditions,  assuming  isentrcpic  expansion  to  sonic 
velocity 

Set _ Ipls 

aw  denotes  adiabatic  wall  conditions 

f  "  flight  conditions 

0  "  reference  value 

1  "  conditions  upstream  of  the  grid 

w  "  conditions  of  the  screen's  wires 
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SECTION  ! 

IN': .,00'JCTI <  ;< 

if  is  well  known  that  by  the  aerodynamic  heating  of  feet ,  vehicles 
flying  at  supersonic,  speeds  arc  heated  to  suc!>  an  extent  that  the 
skin  temperatures  and  the  temperatures  inside  tin*  craft  require  the 
special  attention  of  the  designer.  In  addition,  a  growing  interest  in 
aircraft  vehicle  recovery  operations  has  directed  attention  to  a  variety 
of  deceleration  devices.  In  many  cases  a  sizeable  portion  of  the  initial 
kir -tic  energy  of  the  primary  or  forebody  is  converted  into  an  increase 
in  internal  energy  of  the  retardation  device  and  the  aerodynamic  heat 
transfer  of  the  deceleration  device  becomes  a  design  factor. 

Aerodynamic  decelerators  must  possess  a  large  drag-to-weight  ratio 
and  normally  must  be  capable  of  being  stored  in  a  small  volume.  There¬ 
fore,  most  drag-producing  devices  arc  relatively  thin  and  flexible. 

Towed  decelerators  such  as  ribbon  or  hyperflo  parachutes ,  balloons,  bal- 
lutes,  etc.,  have  flexible  canopies.  Other  proposals  incorporate  strut- 
tural  members  (cones  or  paragliders)  or  hybrid  systems  such  as  the  AVCO 
drag  brake  (1)*  and  the  NASA  rocornet  (2). 

Impermeable  retardation  devices  are  aerody namicaily  unstable  when 
utilized  in  a  supersonic  flow.  As  a  result,  nearly  all  of  the  retarda¬ 
tion  devices  discussed  above  employ  extended  porous  sur i aces  in  con¬ 
junction  witii  venting  techniques  (see  Figure  1).  Most  riven;  1  a.  <  ?  tex¬ 
tile  fibers  either  melt  oi  are  seriously  degraded  by  exposure  to  tempo ra- 

*Numbers  enclosed  by  parentheses  refer  to  references. 
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cures  above  500 °F.  It  is  desirable  to  replace  these  fibers  with  materials 
of  higher  thermal  durability  such  as  fine  metallic  wires  in  which  opera¬ 
tional  temperatures  are  close  to  200G°F. 

The  flow  field  associated  with  a  parachute  results  from  the  inter¬ 
action  of  the  large-scale  phenomena  enveloping  the  complete  body  and 
the  small-scale  processes  associated  with  the  local  flow  through  a  single 
porous  element  in  the  mesh  material  itself.  TK<'  large-scale  phenomena 
(parachute  shape  and  size,  velocity,  altite.de),  which  vary  with  the  par¬ 
ticular  design,  are  generally  involved  with  the  aerodynamics  of  the  prob¬ 
lem  and  contribute  end  or  boundary  conditions  (pressure  ratio,  Reynolds 
number)  to  che  small-scale  phenomena. 

Wnen  opened  at  a  flight  Mach  number  exceeding  unity,  a  parachute 
produces  a  bow  shock  immediately  forward  of  the  canopy.  The  flow  im¬ 
pinging  on  the  parachute  material  then  is  subsonic  and  normal  to  the  sur¬ 
face  except  near  the  periphery  of  the  caonp>  (3).  Hence,  for  most  of  the 
surface  area  of  the  parachute  an  appropriate  model  is  a  woven  mesh  with 
the  impinging  flow  normal  to  it.  The  flow  is  isent Topically  expanded  to 
sonic  onditlons  at  the  minimum  area  of  the  mesh,  based  on  the  open  poros¬ 
ity  of  the  canopy  roof,  and  t  lie  flow  then  expands  from  the  chocked 
condition  into  supersonic  joes  bordered  by  fr*>e  shear  layers  downstream 
(.small-scale  phenomena) .  Immediately  behind  the  individual  fibers  a 
small  region  of  recirculating,  subsonic,  vortex  flow  exists. 

Heat  will  be  transferred  from  the  hot  air  to  the  screen.  The  heat: 
transfer  coefficient  for  a  similar  configuration  is  conventional ly  jx- 


pressed  by  the  following  dimensionless  relationship: 


Nu  »  f(Re,  Pr,  Ma^,  Po)  (1) 

Since  a  sonic  throat  condition  occurs  for  supercritical  pressure 
ratios,  it  appears  of  advantage  to  use  a  sonic  Reynolds  number  defined  as 


Re* 

d 


p*a*d 

V* 


(2) 


where  the  voscosity  y*  is  independent,  of  the  pressure  p*  and  a  known 
function  of  the  temperature  T*. 

The  Prandtl  number  Pr,  whose  value  for  air  is  approximately  0.7 
under  normal  atmospheric  conditions,  deviates  up  to  +  10  percent  from 
this  value  to  temperatures  of  approximately  6000°R,  from  whereon  the  de¬ 
viations  become  larger  (4).  A  temperature  of  6000°R  corresponds  to  the 
stagnation  conditions  at  Ma  =  10  and  50,000  feet  altitude.  At  altitudes 
up  to  400,000  feet  the  stagnation  temperatures  for  Ma  «  10  will  be  below 
6000°R.  The  deviation  of  the  Prandtl  number  from  the  value  0.7  will  be 
neglected  in  this  report. 

The  flight  Mach  number  Ma^  has  only  an  indirect  influence  upon  the 
heat  transfer  to  the  canopy  because,  in  conjunction  with  the  flight 
altitude,  it  will  fix  the  stagnation  pressure  and  temperature  downstream 
of  the  shock.  The.  number  appearing  in  the  relation  for  the  Nusaelt 

number  would  have  to  be  one  Indicating  the  flow  conditions  downstream 
from  the  wire  grid  because  the  upstream  Mach  number  gives  n  >  information 
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on  downstream  conditions  for  a  supercritical  pressure  ratio.  It  is  then 
more  straightforward  to  use  the  pressure  ratio  F  /Pjown  in  place  of  a 
Mach  number. 

•k 

The  porosity  Fo  is  used  as  a  parameter  describin'*,  the  geometry  of 
the  wire  screen. 

The  purpose  of  the  experiments  described  in  this  report  is  to  es¬ 
tablish  the  relationship  between  the  averse*,  Nuaselt  number,  the  Reynolds 

number  Re*  ,  the  pres,  ure  ratio  P  /P.  ,  u»  i  the  porosity  Po  for  a 

d  up  down 

fluid  with  a  Prandtl  number  equal  to  0.7. 


» 
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★ 

Porosity  is 
total  area  of  the 


defined  geometrically  as 
forward- facing  surface  of 


the  ratio  of  open  area  to 
the  retardation  device. 


....  * 


r> 


SECTION  II 


EXPERIMENTAL  FACILITY 

A,  LOW  DENSITY  WIND  TUNNEL 

Figure  2  shows  the  experimental  facility  used.  A  small  subsonic 
wind  tunnel  was  constructed  to  simulate  the  flow  rates  and  reduced 
pressures  characteristic  of  the  operation  of  supersonic  parachutes.  A 
small  duct  2  inches  in  diameter,  in  which  the  flowmeter  was  located,  was 
expanded  into  an  8-inch  O.D.  stilling  chamber  prior  to  flow  through  the 
test  section.  A  mechanical  vacuum  pump-blower  combination,  with  two 
interchangeable  blowers,  provided  the  necessary  flow. 

The  open  system  is  shown  schematically  in  Figure  3.  The  system 
draws  air  from  the  laboratory  room.  It  composed  of  five  components : 
air  dryer,  flow  metering  device,  upstream  and  downstream  valves,  test 
chamber,  and  vacuum  pump  combination,  each  described  in  detail  below. 

1 .  Air  Dryer 

The  air  dryer  is  an  open-ended  cylinder  filled  with  calcium  sul¬ 
fide  and  one  layer  of  glass  wool.  The  dryer  serves  to  dehumidify  the 
incoming  air  below  the  limits  specified  for  proper  pump  operation.  It 
also  provided  a  means  of  eliminating  any  condensation  shocks  that  might 
occur.  The  calcium  sulfide  was  changed  daily  during  running  periods 
to  ensure  optimum  drying. 

2 .  Flow  Metering  Device 

For  the  range  of  Reynolds  numbers  anticipated,  three  different 
metering  devices  were  selected.  Two  sharp-edged  flat-plate  orifices. 


o 


one  1/2  inch  O.D.,  the  other  1/8  inch  O.D.,  and  a  coiled  capillary  tube 
were,  employed.  The  mass  flow  through  the  system  is  found  by  measuring 
the  pressure  drop  across  each  device  when  in  use. 

The  upstream  and  downstream  flanges  were  constructed  so  as  to  per¬ 
mit  interchangeability.  The  corresponding  ranges  for  each  device  are 
listed  below: 


capillary  0.1  <  Re*  <  8 

1/8"  orifice  2  <  Re*  <  100 

1/2"  orifice  80  <  Re*  <  2000 

All  three  flow  metering  devices  were  individually  calibrated.  The 
orifice  plates,  each  constructed  following  A.S.M.E.  specifications,  were 
calibrated  with  water  using  the  weigh-tank  technique.  Calibration  ac¬ 
curacy  was  found  to  be  within  1%.  The  coiled  capillary  tube  was  calibrated 
in  place  using  a  volumetric  flow  apparatus  designed  by  S.  McComas  (5). 

This  calibration  likewise  produced  an  accuracy  of  within  IX. 

The  manometer  board  was  constructed  with  three  U-tubes.  The.  first 
was  filled  with  oil,  giving  the  pressure  drop  between  the  entrance  of 
the  dryer  and  the  flow  metering  device.  The  other  two  were  filled  with 
oil  and  mercury  respectively  to  measure  the  pressure  drop  across  the 
flowmeter.  The  feed  lines  were  fitted  with  valves  to  permit  readings  to 
be  made  with  either  oil  or  mercury.  This  allowed  much  larger  ranges  of 
Reynolds  numbers  to  be  covered  with  each  device. 


3 .  Upstream  and  Downstream  T'alves 

Two  manually  adjusted  Temescal  hlgh-vacuum  gate  valves  were  in¬ 
stalled  in  the  system.  Poth  have  straight  through  connections  designed 
for  minimum  resistance  to  flow  and  low  outgassing  characteristics.  The 
upstream  valve,  two  inches  in  inner  diameter,  was  in*, called  downstream 
of  the  flow  metering  device.  This  valve  was  used  to  throttle  the  flow 
end  provided  the  main  adjustment  for  upstream  pressure.  The  downstream 
valve  could  be  used  to  change  the  pressure  ratio. 

4,  Test  Chamber 

The  test  chamber  is  shown  in  an  exploded  view  in  Figure  4.  The 
chamber  halves  were  cut  from  an  8-inch  O.D.  brass  cylinder,  fitted  with 
flanges  on  each  end.  Each  end  was  then  bolted  to  the  upstream  and 
downstream  valves  respectively .  In  order  to  achieve  uniform  flow,  a 
porous  plate  and  fine  mesh  screen  were  encased  in  the  upstream  end.  The 
central  pai  of  the  porous  plate  was  blanked  off  to  further  generate  uni¬ 
formity  of  flow.  A  mesh  screen  (60  wires  per  inch)  was  inserted  1-1/2 
inches  downstream  to  further  break  up  the  flow  and  reduce  the  scale  of 
turbulence  in  the  chamber. 

The  downstream  half  of  the  chamber  was  fitted  with  wooden  exit 
blocks  to  provide  a  straight  expansion  angle  of  !>°  Leading  to  the  pump. 
This  prevented  an  abrupt  expansion  Just  aft  of  the  test  grid  and  provided 
a  closer  approximation  to  parachute  performance. 

An  aluminum  plat"  with  a  square  cavity  in  the  middle  was  located 
between  the  two  chamber  halves.  The  test  section  was  mounted  in  this 
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TEST  CHAMBER  ASSEMBLY  DETAIL 


cavity  (see  Figure  5).  The  nozzle  block  was  immediately  upstream  of  the 
aluminum  plate  and  attached  to  it  (Figure.  6).  The  nozzle  block  was  ma¬ 
chined  from  a  linen  phenolic  rod,  sanded  and  hand-polished  in  order  to 
keep  the  roughness  as  low  as  possible.  The  primary  concern  in  designing 
the  nozzle  blocks  was  to  minimize  the  boundary  layer  thickness  at  the 
edge  while  preventing  separation.  The  method  of  computing  the  nozzle 
profile  was  to  compute  the  family  of  streamlines  for  the  step  geometry. 
To  calculate  the  boundary  layer  thickness  and  to  get  some  idea  if  sepa¬ 
ration  occurs,  an  integral  solution  to  the  boundary  layer  equations 
was  chosen.  The  method  of  Timman,  given  in  (6),  was  used,  as  it  allowed 
for  the  large  pressure  variations  as  appear  in  a  nozzle.  The  optimum 
shape  produced  a  boundary  layer  thickness  of 

, _  piM 

6(inches)  ■  0.042//Re  Re  - - 

1 

at  the  nozzle  block  exit,  with  p^,  U^,  and  also  evaluated  at  the 
nozzle  block  exit. 

A  total  pressure  probe  was  installed  upstream  of  the  test  grid  in 
the  stilling  chamber,  arid  a  static  pressure  probe  was  inserted  into  the 
wall  of  the  exit  block,  flush  with  the  wall.  The  upstream  operating 
pressure  and  temperature  determined  the  soric  Reynolds  number.  Re*.  In 
addition,  studies  could  be  made  of  the  pressure  drop  across  the  test 
grid  with  the  two  pressure  probes. 


igure  5.  nluffilaum  Plate  with  Test  Section  -md  Housing  mocks 


Figure  6.  Aluminum  Plate  with  Nozzle 


After  sealing  all  external  joints  of  the.  test  chamber  with  Glypta.l, 
a  moderately  hi  t- vacuum  sealant,  the  system  was  leak-tested  and  found 
capable  of  maintaining  a  vacuum  of  less  than  0.1  tarn  of  mercury.  This 
pressure  level  is  one-twentieth  of  that  corresponding  to  the  lowest 
Reynolds  '.umber. 

^ •  Vacuum-Pump  Combination 

Two  vacuum-pump  combinations  were  selected  in  order  to  study  the 

influence  of  the  pressure  ratio  on  the  heat  transfer  coefficient.  The 

first  was  the  Dresser  SSS  120  150  matched  blower  and  pump.  Its  intake 

-2 

operating  range  is  from  10  mm  of  mercury  to  atmospheric  pressure.  The 
pump-blower  combination  maintains  a  constant  pumping  speed  of  110  cu.  ft. 
per  minute  for  inlet  pressures  iuwn  to  1/2  mm  of  mercury,  below  which  the 
pumping  speed  drops  appreciably.  This  combination  is  the  one  which  has 
been  used  mainly  throughout  these  experiments. 

The  second  vacuum-pump  combination  selected  w,:s  the  Dresser  RS-260 
DS-150  matched  blower  and  putnp.  The  pumping  speed  reaches  a  maximum 
of  240  cu.  ft.  per  minute  at  an  inlet  pressure  of  0.5  nun  of  mercury  and 
decreases  to  170  cu.  ft.  per  minute  at  an  inlet  pressure  of  100  mm  of 
mercury;  when  the  inlet,  pressure  decreases  below  0.1  mm  of  mercury,  the 
pumping  speed  diminishes  very  fast  to  reach  the  value  C  at  an  inlet  pres¬ 
sure  of  i  x  1U  arm  of  mercury.  A  Flexishaft  coupling  was  installed 
just  ahead  of  the  pump  (see  Figure  3)  to  eliminate  vibration  in  the  wind 
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B.  TEST  SECTION 


The  test  grid  assembly  is  pictured  in  Figure  5.  The  housing  blocks 
inserted  in  the  square  cavity  of  the  aluminum  plate  were  modeled  of 
Scotchcast  Brand  Resin  No.  250,  capable  of  withstanding  temperatures  of 
250 °F — well  above  operating  temperatures.  The  contact  surface  between 
plastic  and  metal  was  kept  small  to  minimize  conduction  losses.  The 
test  screen  shown  in  Figure  7  was  sandwiched  in  and  Held  by  the  two  halves 
of  each  copper  buss-bar,  the  under  halves  of  which  were  screwed  to  the 

plastic  housing,  and  the  upper  halves  of  which^were  screwed  to  the  lower 

* 

housing.  Electrical  guard  heaters,  consisting  cf  narrow  strips  of  very 
thin  stainless  steel  (0.002"),  were  provided  for  each  buss-bar.  Heaters 

'  v 

were  isolated  from  the  buss-bar  through  thin  asbestos  paper.  Long  copper 
leads  were  slid  through  the  aluminum  plate  and  screwed  to  the  buss-bars 
and  electric  heaters  respectively.  Iron-cons tantan  thermocouples  were' 
inserted  one  in  each  bass-bar.  These  thermocouples,  together  with  the 
electric  heaters,  were  used  to  control  the  temperature  of  the  buss-bars 
such  that  heat  conduction  losses  could  be  eliminated.  The  flow  cross- 
section  was  3/4  x  3/4  inch..  The  square  configuration  was  chosen  to  elimi¬ 
nate  the  extreme  difficulty  of  producing  uniform  heat  generation  in  a 
circular  geometry. 

The  stainless  steel  test  screens  were  cut  from  woven  wire  screens 
and  silver-soldered  to  the  copper  pieces  shown  in  Figure  7.  The  lateral 
sides  of  tile  wires  were  also  silver-soldered  to  provide  paths  of  uniform 
lonkth  ;«ii,  with  this,  of  the  same  electrical  resistance;  thus,  uniform 
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Figure  /.  One  Test  Screen 


heating  could  be  obtained.  Also,  no  difference  in  the  electric  potential 
should  exist  between  the  wires  at  any  location  where  they  touch. 

Stainless  steel  was  selected  because  of  its  relatively  large  re¬ 
sistivity,  its  high  temperature  resistance  coefficient,  and  its  strength, 
which  allow  smaller  deformations  under  the  influence  of  drag  forces. 

The  porosities  were  selected  from  those  commercially  available. 

Three  screens  have  been  used,  with  the  following  properties: 


Screen  1 

Screen  2 

Screen  3 

Diameter 

0.0065" 

0.007" 

0.0055" 

Wires /inch 

70 

60 

80 

Geometric 

Porosity 

29.8% 

33.9% 

31.4% 

Effective 

Porosity* 

0.329 

0.336 

0.325 

* 

This  is 
studies.  The 

a  different  porosity 
effective  po’oslty  C 

used  by  Heinrich  (7) 
is  defined  as 

for  aerodynamic 

with  U  the  average  velocity  of  the  air  through  the  screen  and  V  a  fic¬ 
titious  velocity,  defined  as 

V  -  2&Z 

P 

where  Ap  is  the  pressure  drop  across  the  screen  and  p  the  density  of  the 
air  downstream  of  the  screen. 
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C.  INSTRUMENTATION 


The.  instrumentation  required  for  the  test  apparatus  is  of  three 
types:  temperature,  pressure,  and  electrical. 

Temperatures  were  sensed  by  30  gauge  iron-constantan  thermocouples 
individually  calibrated  to  an  accuracy  of  +p.l°F.  The  calibration 
points  were  fitted  with  least  square  curves  so  that  the  temperature  could 
be  calculated  on  a  computer.  No  significant  error  is  initiated  by  this 
procedure.  A  list  of  vital  thermocouple  locations  is  given  below: 

1.  Upstream  bulk 

2.  Upper  buss-bar 

3.  Lower  buss-bar 

4.  Manometer  board 

The  upstream  total  temperature  was  sensed  by  a  probe  positioned 
two  inches  upstream  of  the  nozzle  block  and  located  in  the  geometric 
center  of  the  stilling  chamber. 

The  total  upstream  pressure  was  measured  by  an  Alphatron  Type  530 
vacuum  gauge  made  by  NRC  Equipment  Corp.  This  instrument  utilizes  an 
ionization  chamber  with  a  radium  source.  It  has  seven  decade  scales 
ranging  from  1000  to  0.001  mm  of  mercury.  The  accuracy  is  2%  of  full 
scale  reading. 

All  electrical  measurements,  including  the  thermocouple  millivolt 
readings,  were  made  with  a  Dymec  2010A  digital  acquisition  system.  This 
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system  measures  analog  data  including  AC  and  DC  voltages.  Accuracy  is 
O.QEof  full  scale  reading.  The  voltages  measured  were  averaged  over  a 
period  of  one  second.  Current  flow  was  found  by  employing  precision 
shunts,  each  Individually  calibrated  before  use. 

Owing  to  the  inductive  effects  created  by  the  woven  mesh  elements, 

DC  power  was  used  for  all  tests.  The  power  was  suppliea  by  two  six-volt 
storage  batteries  wired  in  parallel  for  the  high  Reynolds  numbers  and 
from  a  specially  wired  six-volt  battery,  giving  three  parallel-connected, 
two-volt  cells,  for  the  low  Reynolds  numbers. 

In  order  to  shield  the  test  screen  from  electrical  noises  a  metal 
screen  was  mounted  far  downstream  of  the  test  mesh.  This  screen,  in 
conjunction  with  the  metal  screen  used  to  increase  the  uniformity  of  the 
flow  and  the  pipe  system,  formed  a  closed,  electrically  isolated  volume. 
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SECTION  III 


EXPERIMENTAL  STUDIES 

HEAT  TRANSFER  MEASUREMENTS 

The  goal  of  this  experiment  was  to  find  the  relation  between  the 
average  heat  transfer  coefficient  h  and  the  sonic  Reynolds  number  Re*. 
The  average  heat  transfer  coefficient  is  given  by 


2 

where  Q  is  the  heat  input,  given  by  the  ohmic  heating  1  R,  with  I 
being  the  intensity  of  the  electrical  current  and  the  R  the  electrical 
resistance  of  the  screen,  nnd  A  the  total  area  of  the  wires  forming  the 
test  screen. 

Originally  the  experimental  apparatus  was  designed  to  work  a  slightly 
different  way.  The  governing  equations  are  given  by  a  heat  balance 
at  the  screen, 

I2R  -  hA(tw  -  taw)  ,  (4) 

and  by  the  well-known  relation 

R  -  R_[l  +  a(t  -  t  )]  (5) 

U  W  Wq 

From  the  variables  appearing  in  the  equation,  I  and  R  are  easy 
to  measure,  the  first  one  by  using  a  calibrated  shunt  and  the  second  one 
with  the  help  of  a  standard  resistance.  Further,  with  a  previous  calibra¬ 
tion  of  the  resistance  R  as  a  function  of  the  temperature,  a  coul  ,  be 
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found  and,  knowing  R  and  a,  f.  and  £  could  be  read  from  the  calibration 

w  aw 

curve.  Once  the  temperature  of  the  screen  is  knowr ,  the  electrical  guard 
heaters  could  be  turned  on  and,  with  the  help  of  a  regulable  resistance, 
the  temperature  of  the  buss-bar  could  be  changed  until  the  thermocouples 
inserted  in  them  showed  a  temperature  equal  to  tw-  In  this  way  heat  ex¬ 
changes  by  conduction  between  the  test  mesh  and  the  buss-bars  could  be 
avoided.  Unfortunately,  the  indicated  reference  value  of  the  electrical 
resistance  of  the  test  screens  sometimes  changed  with  time.  However,  no 
changes  in  the  resistance  coefficient  a  were  observed.  The  absolute 
values  of  the  ch  'ges  of  the  resistance  were  within  the  order  of  magni¬ 
tude  of  the  ch  ngea  due  to  a  temperature  variation  between  0  and  5o°F. 

After  checking  contact  resistances,  electrical  noises,  and  ground 
problems,  no  reason  has  been  found  to  explain  the  variations  of  the 
reference  value  of  the  resistance.  Because  of  the  inability  to  establ.it  1 
a  relation  between  the  electrical  resistance  of  the  test  screen  and  its 
temperature,  the  experiment  could  not  be  carried  out  this  way. 

The  observation  that  the  changes  of  the  reference  value  of  the 
resistance  Rq  took  place,  not  continuously,  but  rather  stepwise  at 
a  certain  moment,  permitted  us  the  development  of  a  different  way  to 
perform  the  experiments.  This  method  will  be  explained  in  the  following 
pe ragraphs . 

It  can  be  shown  that  in  the  case  that  the  reference  value  of  the 
resistance  R  does  not  change  during  a  running  period,  the  resistance  K 
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One  disadvantage  of  this  method  is  that  it  does  not  allow  the 

use  of  the  guard  heaters,  so  that  conduction  losses  could  not  be  controlled. 

Therefore^  an  estimate  of  these  losses  was  necessary  (see  Appendix  II). 

The  sonic  Reynolds  number  Re*  was  changed  by  varying  the  upstream 

pressure  (the  air  temperature  in  the  room  was  almost  constant;  thus  the 

sonic  Reynolds  number  wa  a  linear  function  of  the  upstream  pressure  only) . 

»he  variation  of  the  upstream  pressure  P  was  achieved  by  manipulating 

the  upstream  valve.  When  the  pump-blower  combination  RSS-120-150  was 

used,  th  ’  pressure  ratio  P  /P .  could  be  ma  intained  almost  constant  for 

up  down 

all  th ree  screens — independent  of  che  upstream  pressure  ? — as  shown  in 
Figures  9,  10,  and  11.  The  small  oscillations  shown  may  be  due  to  small 
changes  iu  the  efficiency  of  the  pump-compressor  system  with  the  input 
pressure  and  to  the  small  errors  in  readings  of  the  upstream  and  down¬ 
stream  pressures.  When  the  pump-blower  combination  RS  260  OS- 150  is 
used  for  „.ow  upstream  pressures  as  shown  in  Figure  10,  the  efficiency 
of  the  pump-compressor  system  snows  a  much  larger  sensitivity  to  the  inlet 
pressure,  as  described  in  Section  i  I .  A . . 


SECTION  IV 


DISCUSSION  OF  RESULTS 


The  dimensionless  llusselt  number, 


Ni.' 


* 


(7) 


based  on  the  wire  diameter,  has  been  used  instead  of  the  dimension- 
dependent  heat  transfer  coefficient.  The  results  experimentally  ob¬ 
tained  are  shown  in  Figures  12,  13,  and  14,  where  the  Nusselt  number 
Nu  has  been  plotted  as  a  function  of  the  sonic  Reynolds  number  Re*.  In 
Figure  13  the  values  obtained  by  using  the  higher  pressure  ratio  are 
also  included.  Although  the  influence  of  the  pressure  ratio  ^p/^jcwro 
appears  to  be  small,  a  tendency  towards  decreasing  the  dimensionless 
heat  transfer  coefficient  with  increasing  pressure  ratio  can  be  observed. 
This  phenomenon  seems  to  be  logical  since  a  decrease  of  the  pressure 
ratio,  while  holding  the  upstream  pressure  constant,  indicates  a  lower 
downstream  static  pleasure  which  should  affect  the  local  heat  transfer 
coefficient  in  the  manner  noted. 

No  effects  of  the  porosity  Fo  on  the  average  Nusselt  number- — within 
the  relatively  narrow  range  of  poroeities  investigated — can  be  observed 
irt  Figures  12  through  14. 

It  is  of  interest  to  compare  the  measured  Nusselt  numbers  for  a 
scr-e o  with  Nusselt  numbers  obtained  for  f  low  across  a  cylinder.  Figure 
15  t  the  result  of  such  a  comparison.  The  lines  in  the  figure  represent 
average  Nusselt  numbers  as  a  function  of  upstrectn  Reynolds  number,  with 
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Figure  15. 


Comparison  of  Present  Results 
with  Results  for  a  Cylinder 


the  upstream  Mach  nuitber  as  a  parameter.  These  Nusselt  numbers  have 

been  taken  from  (8).  The  symbols  in  the  figure  refer  to  the  results 

of  oar  measurements „  For  this  purpose  the  Reynolds  number  Re^,  based 

on  flow  conditions  at  the  nozzle  block  exit  and  on  the  wire  diameter, 

has  been  obtained  from  the  sonic  Reynolds  number  Re*.  The  Mr eh  number 

* 

Maj  warn  also  calculated  at  this  location.  The  Knuds en  number  calculated 
from  the  relation 


✓foe” 


(8) 


is  also  indicated  in  Figure  15  as  abscissa. 

One  should  be  careful  in  drawing  conclusions  from  a  quantitative 
comparison  of  the  curves  with  the  lines  in  this  figure  because  of  the 
difference  in  flow  condition  on  a  single  cylinder  and  on  a  screen.  For 
the  same  upstream  flow  condition,  for  instance,  the  maximum  velocity  out¬ 
side  the  boundary  layer  on  the  cylinder  is  larger  for  the  screen  than  for 
the  single  cylinder.  Nevertheless  a  qualitative  comparison  certainly 
can  be  made  and  the  figure  indicates  that  rarefaction  effects  become  im¬ 
portant  for  Reynolds  numbers  lower  than  a  value  of  order  100.  This  rare¬ 
faction  effect  alr.o  explains  why  tire  Nusselt  numbers  continuously  de¬ 
crease  at  about  the  same  rate  with  decreasing  Reynolds  number,  whereas  in 
incompressible  flow  (Ma  -  0)  the  decrease  of  the  Nusselt  number  diminishes 


For  plotting  purposes  in  Figure  15,  the  Mach  number  Ma^  has  been 
considered  constant  and  equal  to  0.185. 
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as  Che  Reynolds  number  approaches  the  value  1.  The  figure  shows  that 
actually  the  agreement  between  the  Nusselt  numbers  measured  on  the  screens 
will  values  for  a  single  cylinder  is  surprisingly  goods  Only  at  the 
highest  Reynolds  numbers  investigated  (beyond  a  value  Re^  -  250)  a 
sudden  Increase  in  Nusselt  number  Is  indicated  which  is  not  evident 
on  the  curves  for  the  single  cylinder,. 

1/2  1/2 

In  Figure  16  has  been  plotted  the  ratio  Hu/ (Re*)  versus  (Re*) 

1/2 

Here  the  value  Nu/(Ue*)  seems  to  be  constant  for  Re*  larger  than  400. 

In  this  figure  are  also  included  the  integrated  results  reported  in  (9). 

In  this  study  the  local  measu  ements  were  made  on  an  enlarged  woven 
mesh  in  which  the  individual  fibers  were  one  Inch  In  diameter.  The 
porosity  and  pressure  ratio  of  the  data  in  (9)  closely  match  the  present 
studies,  while  the  sonic  Reynolds  number  was  consider anly  nxgher.  The 
agree. ..«_nt  between  the  integrated  local  values  of  (9)  and  the  present 
average  values  is  very  good. 

C.  Scott  derived  the  following  expression  describing  the  Nusselt 
number  at  the  stagnation  point  of  the  wires  In  a  screen  from  stagnation 
point  boundary  layer  analysis  (10)  ; 


In  Figure  11  this  relation  is  compared  with  the  measurements  on  the  three 
screens.  For  sonic  Reynolds  numbers  Re*  above  a  value  of  400,  the  average 


i  r> 


Thcaratlcil  Pr*^icfcion 


Nusselt  number  is  bv  a  factor  of  2  lower  than  he  one  given  in  Eq.  9. 
It  can  therefore  be  recommended  to  use  the  equation 


(L+A) 


Re* 


0.5 


(10) 


tor  the  aiculati  ..  of  average  Nusselt  numbers  of  screens  with  the  porosity 
range  .investigated  in  our  experiments.  For  sonic  Reynolds  numbers  below  400 
the  averagt  Nusselt  numbers  start  dropping  off  with  decreasing  Reynolds 
number  as  a  consequence  of  the  rarefaction  effects. 

Figure  17  also  contains  the  integrated  values  of  (9).  The  agreement 
with  the  present  values  is  also  quite  heartening. 


SECTION  V 


CONCLUSIONS 

Heat  transfer  experiments  have  been  performed  on  three  woven  wire 
mesh  screens  simulati'  ;  porous  parachute  cloth  with  porosities  of  29.8, 
31.4,  and  33.9  *-■  cecal.  Average  Nusselt  numbers  were  measured  as  a 
function  of  sonic  Reynolds  number  at  a  ratio  2.5  of  the  upstream  total 
to  the  downstream  static  pressure.  The  results  are  presented  in  Figures 
12  through  16.  They  indicate  that  rarefaction  effects  start  to  influence 
heat  transfer  when  the  Reynolds  number  Re*  based  on  so  conditions 
decreases  beyond  the  value  of  approximately  400  or  for  a  corresponding 
value  of  the  Knudsen  number.  Within  the  range  of  porosities  which  were 
investigated  and  the  accuracy  obtained,  no  effect  of  porosity  on  the 
average  heat  transfer  coefficient  could  be  established.  The  method  of 
the  ratio  of  the  upstream  total  oressure  to  the  downstream  static  pressure, 
was  checked  on  one  of  the  screens  and  was  found  to  be  small  for  values 
between  2.5  and  4,  The  Nusselt  numbers  for  screens  are  shown  to  have  a 
behavior  quite  similar  to  the  N"sselt  number  for  a  single  cylinder  in 
cross  flow.  The  relation  Nu(L  +  i))/Re*  =  0.5  appears  to  be  an  acceptable 
repre.si  nt  at  ion  of  tin  average  Nusselt.  numbers  on  »<•  reer.o  lor  sonic  Reynolds 
numhei  s  above  400  tor  which  r.»i  e  faction  effects  are 


not  present. 


APPENDIX  I 


DATA  REDUCTION  TECHNIQUE 


Section  III  shows  that  the  governing  equations  are 


I  R  *  hA(t  -  t  ) 
w  aw 


and 

R  -  R0U  +  <.<£„  -  t^)] 

When  the  wind  tunnel  is  running  and  1*0,  then  t  *  t  by  the 

aw  w0 

definition  of  the  recovery  temperature,  so  the  equations  can  be  rewritten 


2 

I  R 


hAt 


(11) 


R-  RC 

V 


At 


(12) 


and,  dividing  Equation  11  by  Equation  12, 


I2RR  a 
_ 0 

A(R  -  R  ) 


(13) 


which  Is  the  formula  used  to  calculate  the  heat  transfer  coefficient. 

In  order  to  show  that  the  resistance  R  is  a  linear  function  of  the  square 
of  the  intensity  of  the  current.  Equation  IS  is  written 


2  hA(R  "  V 


Ry  <1  +  lit) 


(14) 
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Now  in  our  experiments  the  product  aAt  <<  1  (in  general  for  the 
highest  values  of  I,  aAt.  <  0.03)  so  that  this  term  can  he  and 

Equation  14  gives 


I2  *  —  (R  -  R0) 
R0 


which  shows  the  desired  relation  between  1  and  R. 


(15) 
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APPENDIX  II 


ESTIMATE  OF  CONDUCTION  LOSSES 

Because  conduction  losses  could  not  be  eliminated,  an  estimate  of 
them  seems  to  be  important. 

The  model  used  for  the  estimate  of  the  heat  conduction  losses  was 
a  wire  with  uniform  heat  sources  stretched  between  two  walls,  of  prescribed 
temperature,  with  the  intensity  of  the  heat  source  given  by  the  ohmic 
heating.  The  solution  is  given  in  (11).  The  wire  is  cooled  by  con¬ 
vection.  Its  length  is  equal  to  the  actual  length  of  the  wire  measured 
from  one  buss-bar  to  the  other  across  the  soldered  connection  (Fig.  7). 

The  average  heat  transfer  coefficient  measured  in  the  experiments  was 
used  in  the  analysis.  In  Figure  18  are  presented  the  results  of  the  heat 
conduction  correction  for  some  key  points. 

An  estimate  of  heat  losses  by  radiation  was  also  made  assuming 
the  shape  factor  for  this  exchange  of  radiant  energy  to  be  one.  The 
heat  transferred  by  radiation  in  this  case  will  be  larger  than  that  cor¬ 
responding  to  the  actual  situation.  As  a  result  we  saw  that  the  heat 
losses  by  radiation  are  of  one  order  of  magnitude  smaller  than  the  heat 
conduction  losses,  so  they  could  be  completely  neglected. 


EXPERIMENTAL  RESULTS 
RES  ULTS  CORRECTED  FOR 
CONDUCTION  LOSSES 
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